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X ■ ABSTRACT 

We calculate the energy-momentum distribution of the charmed quarks produced 
in neutrino reactions on protons, quantifying the importance of mass and current non- 
conservation effects. We study the strange and charm distributions probed in neutrino 
interactions in the presently accessible kinematical region. Some ambiguities inherent 
to the extraction of the parton densities from dimuon data are pointed out. 



Charm excitation in neutrino scattering is a primary source of information on the 
strangeness in nucleons. At moderate values of Q 2 (the momentum transfer squared 
in deep inelastic scattering), subtle threshold and current non-conservation effects are 
at work (I], |2|, [3], |J, which lead to the non-universality of the charm contribution 
to the structure functions probed by muons and by neutrinos. The driving term of 
the neutrino-excitation of charm is the W-gluon fusion process W + g — > cs (Fig. 1), 
which consists of two partonic subprocesses: excitation of charm on the strange sea 
and excitation of anti-strangeness on the anti-charmed sea (we shall consider only the 
Cabibbo unsuppressed reactions). Both subprocesses are of course an integral part of 
the gauge-invariant QCD cross section and gauge invariance requires both amplitudes 
(Figs, la, lb) to be taken into account. In principle, the two subprocesses can be 
separated according to which parton, c or s, is produced in the W^-boson emisphere in 
the Breit reference frame (other conventions are possible, as well). Such a separation 
becomes well defined at Q 2 ^> (/i + m) 2 , which is not yet the case with the present 
experiments (hereafter, m and [i will denote the mass of the strange and of the charmed 
quark, respectively). 

From an experimental point of view, the extraction of what is called the strange 
structure function of the nucleon is a multi-stage process (see || and the references 
therein). The signature of the excitation of the charmed quark is the semileptonic decay 
c — > s/i + z/, which leads to the opposite-sign muon pairs in the final state (here /i~ 

comes from the primary weak interaction vertex). The main source of background are 
the semileptonic decays of pions and kaons. This forces one to put stringent lower cuts 
on the energy of fi + . Because of these cuts the measured cross section is sensitive to the 
energy of the produced charmed quarks. In the limit where the masses of the quarks 
can be neglected, the primary produced quark carries ~ 100% of the W's energy, and 
then fragments with some fragmentation function (usually supposed to be independent 
of Q 2 ) into the charmed particles, the semileptonic decays of which produce the de- 
sired signature of the charm production. At moderate values of Q 2 , the mass effects 
deeply modify this scenario. In this paper we shall show that, at moderate Q 2 , the 
primary charmed quarks have quite a broad energy distribution. Furthermore, this 
energy distribution depends on Q 2 , and understanding this dependence is important 
for unfolding the acceptance effects and for the determination of the primary charm 
production rate. 

We shall focus on the region of small values of the Bjorken variable x = Q 2 /2m p v 
and on the lowest-order QCD subprocesses, which simplifies the analysis. However, it 
will be evident that our main conclusions have a broader applicability. 

We have to calculate the W- and 7*— absorption cross section daT,L/ 'dzcPk for 
the transverse (T) and longitudinal (L) W and 7*. We denote by k the transverse 
momentum of the charmed quark, and by 



z = 




(1) 



2 



the fraction of the light-cone momentum of the W (7*) carried by the charmed quark. 
Here we concentrate on the (z, k) distribution of the produced quarks, i.e. of the jets 
generated by these quarks, leaving aside the issue of quark fragmentation into charmed 
hadrons. 

The analysis of the energy and momentum distribution in the case of massless 
quarks and/or asymptotically large Q 2 was first performed by Altarelli and Martinelli 
0. More recently Ellis and Nason [|7| gave a very detailed treatment of the real 
photoproduction of charm. One important finding was that the higher-order QCD 
corrections only slightly change the energy and transverse momentum distribution of 
the produced quarks. 

In this paper we shall be concerned with the region of moderate Q 2 relevant to 
neutrino experiments, where the strong unequality of masses of the charmed and 
strange quarks manifests its effects. Two important points are the onset of the par- 
ton model interpretation of heavy flavors when passing over the threshold value of 
Q 2 ~ 4/x 2 , (m+/i) 2 , and the transition from the dominance of the longitudinal cross sec- 
tion and the strong breaking of the Callan-Gross relation at small Q 2 to the dominance 
of the transverse cross section and the (approximate) restoration of the Callan-Gross 
relation at large Q 2 . 

To the lowest order in QCD, the longitudinal and transverse absorption cross sec- 
tions in the charged current case are given by a generalization of the formulas provided 
in m - see eqs. (11) and (12) - and read [a ew is the appropriate electroweak coupling 
constant for the specific process considered) 
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where the gluon-gluon-nucleon vertex function V(k) is related to the charge form factor 
of the proton G em (q 2 ) through V(k) = 1 — G em {?>n 2 ), 

e 2 = z(l - z)Q 2 + zm 2 + (1 - z)fx 2 , (3) 

and the functions ^(fc^ fc 2 ) are explicitly given by 

Nrik^kz) = [z 2 + {l-z) 2 }{g 2 v + g 2 A )k l -k 2 

+gl[zm + (1 - z)fi} 2 + g 2 A [zm - (1 - z)fi] 2 , (4) 



N L (ki,k 2 ) = — ([gy{m - fi) 2 + g 2 A (m + ■ k 2 

+g 2 v {2Q 2 z{l -z) + (m- n)[zm - (1 - z)fx}} 2 

+ g 2 A {2Q 2 z(l-z) + (m + ij)[zm + (l-z)fi}} 2 ) . (5) 
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In muon scattering gy = (the charge of the quark involved in units of the electron 
charge), = and fi = m. In the charged current (CC) neutrino interactions gA = 
—gy = — 1 and m and fi stand for the strange and the charm quark masses. In the 
neutral current (NC) neutrino interactions fi = m and the corresponding vector and 
axial coupling are given by the Standard Model ||. The strong coupling as{k 2 ) in 
front of the differential cross section enters at the virtuality of the quark k 2 = e 2 + k 2 . 

The ^-distributions of the charmed quark can be calculated by integrating eq. (Q) 
over the transverse momentum. 

Let us start with the ^-distribution for <jt- Here Nt is a mild function of z, and 
the z dependence comes from e 2 . To a crude approximation, neglecting the scaling 
violations, we find 

dor z 2 + (l-z) 2 = z 2 + (l-z) 2 

dz ~ e 2 z(l - z)Q 2 + zm 2 + (1 - z)fi 2 ' [ ' 

Consider first the muon scattering, fi — m. The distribution (|6|) is nearly flat at 
Q 2 < 4/i 2 ~ 10 GeV 2 /c 2 . Only at asymptotically high Q 2 > 4/r 2 , it develops the 
parton model peaks at z — > and z — > 1 (see 0), so that the charmed (anti)quark 
carries a fraction z ~ l—fi 2 /Q 2 of the photon's light-cone momentum. The result of the 
exact calculation of (l/a^dax/dz for the cc excitation in electromagnetic scattering is 
shown in Fig. 2a. The situation is similar for the NC neutrino scattering (Fig. 2b). 

In the CC neutrino scattering the charmed quark is much heavier than the strange 
quark, fi ^> m. At large Q 2 the z distribution (shown in Fig. 2c) again develops peaks at 
z — * 1 and z — » 0. However, the Q 2 evolution of the forward, dax/dz ~ (1 — z— m 2 /Q 2 ), 
and the backward, dar/dz ~ (1 — z — fi 2 /Q 2 ), peaks is quite different: the forward 
peak appears sooner. 

Now, the CC excitation of charm and (anti) strangeness are inseparable. In the 
parton model language, at very large Q 2 , the peak at z — > 1 can be identified with 
the excitation of the charm on the strange sea, while the peak at z — * describes 
the excitation of s on the c distribution in the nucleon. Thus, only asymptotically the 
separation of the two subprocesses depicted in Fig. 1 is well defined. As we have already 
mentioned, in the extraction of the dimuon data, a cut on z is implicitly made, which 
produces an acceptance-dependent separation of the strange and charm contributions 
to the W-absorption cross sections. However, at moderate Q 2 the sensitivity to the 
cutoff z c is rather strong and the parton model reinterpretation of the cs excitation 
in terms of the two partonic subprocesses is by no means unique and becomes matter 
of convention. In this region of Q 2 an assumption of the form da /dz ~ 5(z — 1), in 
conjunction with a Q 2 -independent fragmentation function, is untenable because it 
would neglect the crucial Q 2 -dependence of the ^-distributions. 

The problem is complicated by the presence of the longitudinal contribution, which 
contribute significantly to the "non-partonic" domain of z around 1/2 (for a discussion 
of the massless case see 0). 



4 



In muon interactions, again neglecting the scaling violations, one finds a broad 
symmetric z-distribution, shown in Fig. 2a 

da L Q 2 z 2 (l-z) 2 
dz ~ 5 4 

which, in proximity of the forward peak, increases with Q 2 . The ratio R = Ohl a T is of 
course rather small, R ~ Q 2 /4fj, 2 . 

By contrast, in neutrino scattering, because of the nonconservation of the weak 
axial current (gA 7^ 0) and of the flavor changing vector current (fi 7^ m), one has 
N L r*j Nt^i 2 /Q 2 and the longitudinal contribution is rather significant at moderate Q 2 
(R ~ 4fi 2 /Q 2 ). In the NC interactions the ^-distribution is symmetric around z — ^ 
(Fig. 2b). In the CC interactions, da^/dz develops two peaks at z — > 1 and z — ► 0, 
and is asymmetric around z — | (Fig. 2c). 

An interesting quantity is -R(z) = <Jl{z) / &t{z) which shows how the weak-current 
non-conservation effects depend on z. In Fig. 3 we compare R(z) at Q 2 = 4 GeV 2 /c 2 for 
cc excitation in muon scattering and in NC neutrino scattering, and for cs excitation in 
CC neutrino scattering. Shown for comparison are the values of R for the ^-integrated 
cross sections. In virtual photoabsorption both R and R(z) are small. In the NC and 
CC neutrino scattering R(z) decreases as z, (1 — z) — > 0. Notice that in the CC case 
R(z) is slightly asymmetric. 

The peculiar Q 2 -dependence of the transverse and longitudinal contributions to 
the ^-distributions manifests itself in a detectable way at the level of the sea parton 
densities. We must say that neither the introduction of the density of partons at 
Q 2 < 4/i 2 , (m + yU 2 ) 2 , nor the comparison of the ar-dominated muoproduction with 
the cTL-dominated neutrino-production do make much sense. Nonetheless, in order to 
clarify the issue of the non-universality, let us proceed with such a comparison. 

The well-defined quantities are the cross sections (Jt,l, which can be converted into 
the structure functions F LjT (x,Q 2 ) = Q 2 a^^l^^ew In terms of the more familiar 
structure functions Fi,F 2 , one has F 1 = F T /2x and F 2 = F T + F L . At very large Q 2 , 
when all quarks can be regarded as massless, the structure functions can be decomposed 
in terms of the parton densities as M (we consider an isoscalar nucleon) 

i \ 4 1 - 

F 2 W {x, Q 2 ) = - [u„ + u M + c M + c M ] + - [dp + d^ + s M + s M ] (8) 

F 2 {u \x, Q 2 ) = u v + u v + 2 c v + d v + d„ + 2 s u (9) 

These formulas can be taken as the operational definition of the parton densities even 
at Q 2 < 4// 2 , (m + fi) 2 . Similarly, one can introduce the parton densities defined 
in terms of Ft(x,Q 2 ), which would have been identical to the ones in Eqs. (fH]) were 
it not for the breaking of the Callan-Gross relation. Notice that we have supplied the 
parton densities by subscripts v and \i. The non-universality of the charm and strange 
densities is quantified by the ratio r u /^ = (c u + s y )/(c M + s M ) (however one has to keep 
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in mind that there is no direct measurement of We have evaluated r u /^ within 
the model of Ref. [g § 0, g, |IDf taking (m + p) 2 = AGeV 2 /c 2 . The result is shown in 
Fig. 4a, where both the x- and the (^-dependence of r v /^ is exhibited. Notice that, 
whereas at large Q 2 the ratio r v / il tends to the asymptotical (and naively assumed) 
value of unity, at small Q 2 (< 10 GeV 2 /c 2 ) it increases with Q 2 before flattening down. 
In the kinematical region probed by the CCFR experiment |J the deviation of r u / fJi from 
unity is as large as ~ 25% and is a possible explanation of the observed discrepancy 
between different determinations of the strange density |f2j . The numerator in r v /^ is 
dominated by the longitudinal component: taking only the transverse contribution, 
the departure of r u /^ from unity would be even stronger. For completeness we present 
in Fig. 4b the ratio ^/(c^ + s^) evaluated at two different values of Q 2 . 

Whereas in the determination of F 2 one deals with the sum c v + s u , which is free 
from ambiguities, being related to the integral of d<jT,L/dz over the whole z-range, the 
experimental analysis of the dimuon data introduces a separation of s u and c u . The 
cs production cross section is thus splitted between the two partonic subprocesses of 
Fig. 1. This amounts to introducing a cutoff z c in the ^-distributions, so that s u and c v 
are defined in terms of the cs production cross section subject to the (arbitrary) cuts 
z > z c and z < z c , respectively. In the treatment of the dimuon data the cutoff z c is 
implicitly posited and the strange density consequently obtained. Were the data taken 
at asymptotically large Q 2 , when the peak at z = 1 is delta-like, the choice of the 
cutoff would be irrelevant. However, the Q 2 region of interest, especially at small x, is 
that of small and moderate values Q 2 < 20 GeV 2 /c 2 , where the choice of z c necessarily 
introduces some arbitrariness. 

Choosing the cutoff z c = 1/2 is equivalent, in the muon interaction and in the 
weak NC interaction, to simply equating the charm and anticharm densities, since the 
^-distributions are symmetric. Also in the CC interaction at asymptotically large Q 2 , 
where the two parton model peaks tend to become delta-like, taking z c = 1/2 gives 
s v — c v — (c v + s v )/2. By contrast, in the CC case at small Q 2 the ^-distributions 
are strongly asymmetric , the backward and the forward peaks are barely visible and 
the splitting of s u and c v cannot be unambiguously done. Hence some care must be 
exerted in analyzing the dimuon data in terms of the strange density. The effect of two 
different choices of the cutoff z c on s u [z c = 1/2 and z c = 0.8, the latter corresponding 
approximately to taking s u = c u = (s u + c u )/2) is displayed in Fig. 5. Our curves are 
presented for Q 2 = lOGeV 2 /c 2 and compared to the dimuon data around the same 
value of Q 2 . An exact comparison with the data would require the knowledge of the 
experimental value of z c , which is lacking. 

In conclusion, the calculation of the energy-momentum distribution of the quarks 
produced in the charm leptoproduction has shown that, in the presently explored 
kinematical region, mass and current non-conservation effects play a decisive role, 
making the extraction of the heavy-quark parton densities from dimuon data subject 
to some (solvable) ambiguity. We have expressed in a quantitative way the deviations 
from the naive parton model expectations and stressed the importance of a correct 
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interpretation of the data coming from neutrino interactions. 
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Figure captions 



Fig. 1 - The excitation of charm in the FT-gluon fusion process. 

Fig. 2 - a) The ^-distributions (l/cr TtL )daT t L/dz of the charmed quark produced in the 
transverse (T) and the longitudinal (L) virtual photoabsorption. The curves are 
as follows: T at lGeV 2 /c 2 , dotdashed; T at 10 GeV 2 /c 2 , solid; L at lGeV 2 /c 2 , 
dotted; L at 10GeV 2 /c 2 , dashed. 

b) Same as a), but for the weak neutral current case (Z°-absorption cross sec- 
tions). 

c) Same as a), but for the charged current case (W-absorption cross sections). 

Fig. 3- The ratio R(z) = (Jl{z)/(Jt{z) at Q 2 = 4GeV 2 /c 2 in muon scattering (solid 
curve), CC neutrino scattering (dashed curve), NC neutrino scattering (dot- 
dashed curve). 

Fig. 4 - a) The ratio (c„ + s„)/(c M + s M ) plotted as a function of x at three different 
values of Q 2 (solid curve, AGeV 2 /c 2 ; dashed curve, lOGeV 2 /c 2 ; dotdashed curve, 
30GeV 2 /c 2 ). 

b) The ratio c M / (c^ + s M ) at two different valus of Q 2 : 10 GeV 2 / c 2 (dashed curve), 
30 GeV 2 /c 2 (dotdashed curve). 

Fig. 5 - The strange density s v probed in neutrino interactions calculated at Q 2 = 10 GeV 2 / c 2 
with two different cutoffs: z c = 0.5, solid curve (T+L), dashed curve (only T); 
z c = 0.8, dotdashed curve (T+L), dotted curve (only T). The data are from 
Ref. [5]. 
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